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THE LINEAR DIFFERENCE EQUATION OF THE FIRST ORDER.* 

By K. p. Williams. 

In this paper we consider the first order linear difference equation 

(A) fix + 1) = E(x)f(x), 

where E(x) is an entire function. It will be seen that the theory for such 
an equation gives at once the theory for the equation where the coefficient 
of /(x) is a meromorphic function. 

The equation (A) was first completely solved by Hurwitz,t who showed 
the existence of a solution when E{x) is any entire function. He did not, 
however, consider this equation primarily, but rather the equation 

(B) fix + 1) - m = <p{x), 

for which he proved that there was an entire solution in case <p(x) is entire,| 
and a meromorphic solution in case <p(x) is meromorphic. The logarithmic 
derivative of the first equation reduces to one of the form (B), where the 
known function is meromorphic. Thus the solution Hurwitz gave of (B) 
included that of our equation. 

Somewhat later Barnes§ gave a direct treatment of equation (A), by 
expressing the function E{x) in its Weierstrass factor fonn. To include all 
cases he found it necessary, however, to know that equation (B) has a 
solution when <p(x) is an entire function. But Hurwitz proved in an elegant 
way that to pass from the case where <p{x) in (B) is entire to the case where 
it is meromorphic, is but a simple step, and this, as remarked above, is to 
know completely the theory for (A). Thus we gain nothing if we appeal 
at all to (jB) in giving a direct solution for the first equation. The solution 
that Barnes obtained is, furthermore, to a large extent undetermined; for 
it contains certain functions of the zeros of E{x), whose explicit forms, he 
states, he was able to determine only in certain simple cases. 

The present discussion also depends fundamentally on the factor form 
of E{x) ; but the method is essentially different from that used by Barnes. 
The solution which is obtained is simple, and is completely determined. 

* Presented to the American Mathematical Society, December 27, 1913. 

t Acta Mathematica (1897), T. XX, pp. 285-312. 

t An interesting discussion of (B) for the case where <p(x) is entire is given by Carmichael, 
Amer. Journal of Math., vol. XXXV, pp. 165-182. 

§ Proceedings of London Math. Soc. (1904), 2d series, vol. 2, pp. 438-469. Barnes considers 
the equation f(x + u) = <p{x)f{x), where w is any real or complex constant. This is no more 
general, and adds only to the complexity of the formulae. 
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130 K. p. WILLIAMS. 

We restrict ourselves entirely to those cases which do not need to be made 
to depend on (5) ; this includes all cases where E{x) is an entire function 
of finite genus, and also a large class of cases where this limitation is not 
fulfilled. It is also thought that the present method has the advantage of 
showing clearly why a simple choice of certain convergence factors can be 
made. 

Let ai, tti, • • • , be the zeros of E(x) ; then we know by Weierstrass's 
factor theorem that we can always write 

1-1 , 



E(x) = e''<^'Il(l--)e"='"^'"^ 

t=l \ di / 



where g(x) is itself some entire function. If we can find an integer k such 
that the series 

« 1 k+l 

X - 

is convergent, we say the function E{x) is of finite order, and we have the 
simpler form 






If further g{x) reduces to a polynomial of degree q, we say that the function 
E{x) is of finite genus r, where r is the greatest of k and q, provided k is the 
smallest positive integer for which the above series is convergent. We shall 
limit ourselves to the cases where g{x) is a polynomial, the order being either 
finite or infinite. 

The resolution of E{x) into factors suggests at once a method for solving 
our equation; for a fundamental property of such an equation is that the 
solution of 

/(x + 1) = Ui(x)U2{x)f(x), 

where Ui{x) and Ui{x) are known, is found by forming the product of the 
solutions of the two equations 

fix + 1) = Ui(x)f(x), and /(x + 1) = Ui(x)f{x). 

In order to be able to write the solutions corresponding to the exponential 
factors in ^(x) we note first that a solution of the simple equation 

fix + 1) = e-"/(x), 
where a is any constant, and n a positive integer, is given immediately by 

fix) = e "•+' , 
where <Pn+ii^) is the nth Bernoulli polynomial. The fundamental prin- 
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ciple mentioned above then shows at once how to solve the equation 
f{x + 1) = e'^^^fix), where v(x) is any polynomial. In the cases we are 
considering we can thus solve easily the equation 

fix + 1) = e"(-'/(x), 

which arises from the extraneous exponential factor in E(x); we indicate 
the solution by e''-''\ g(x) being a polynomial which can be expUcitly ex- 
pressed in tenns of the BernoulU poljmomials and the constants in g(x). 

If now we assume that equation (B) has a solution whenever (p{x) is 
entire, we see that the equation which we have last written has a solution 
even in case g(x) ceases to be a potynomial, but is any undetermined entire 
function. It was precisely at this point that Barnes introduced (B) into 
his consideration of (A). 

We notice finally that a solution of the equation f{x+ 1) = (1 —x/ai)f(x), 
is given by/(a;) = T(x — a,)(— a,)"^, r(x) being the gamma function. We 
thus have solutions corresponding to all the types of factors in Eix) ; and 
from these we build up a solution of (A). 

Let us first assume that the function E{x) is of finite order. The 
equation then takes the form 

fix + 1) = /(x)e«<^' n { 1 - ^ ) e"=' " ^'" , 

being convergent. 
A formal solution is at once seen to be 



°° I 1 
the series 2 ~ 

7^1 \ at 



where 






o r^\ 'P2(^) _i_ ^3(3;) <Pk+i(x) 



1 ■2ai ' 2 • 3a<2 ' ^ kik + l)a<* ' 



while pi, P2, • • • are any constants (or periodic functions of period 1). 

If now we can choose the quantities p, so that the above expression is a 
uniformly convergent product we shall have an analytic solution of our 
equation. A method for doing this will be found by means of the asymp- 
totic form of the gamma function. 

Before we proceed, however, we shall introduce a certain convention 
which will save unnecessary repetition. We have frequent occasion to 
use functions of x and a<, say M(x, Oi), which have the property of being 
less in absolute value than some properly chosen constant, provided x is 
confined to any Umited region, and \ai\ is large enough. Since |a,| in- 
creases indefinitely with i, the last condition will be fulfilled if i is sufficiently 
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great, say if i > /. The value of / depends in general on the region in 
which X can vary. Under these conditions we shall merely say that 
M{x, ai) is bounded, always understanding by that term that x and i satisfy 
the requisite conditions; that is, x is in some limited region, while i is larger 
than some positive integer, depending perhaps on the region in which x 
can vary. Similarly, when we say a function of a,- alone is bounded, we 
shall always mean that it is less in absolute value than some constant if 
|a,| is sufficiently great, that is, if we choose i large enough. Such functions 
as those described will occur in our work upon breaking off convergent 
series, or series which asymptotically represent a function. 

In order to be able to expand T{x — a,), we shall assume for the present 
that all the points Ci, Cj, • • • lie to the left of a certain line parallel to the 
imaginary axis, and are arranged in order of increasing absolute values. We 
can then use Stirling's formula for the gamma function, and we thus obtain 

r(x - a.) = V2ir(x - ai)^^<-ie-^^'+*<'^'«', 
where 

Bi Bi 



^(X - Ci) = 



1 • 2(x - a.) 3 • 4(x - UiY 



(-ly+'Bt Rt(x,a,) 

■^ ^ {2t- l)2<(x - aiY*-^ "•" (x -aO^'+i ' 

the B's being the BemouUi numbers, t any positive integer, and Bt{x, a,) 
a bounded fimction. 

Therefore if i is sufficiently great we have for the general term of our 
formal product 

r(x - o.-)e*"'> ^ V2t( - a.)-°^> ^,,, 
i-a,ypi p, ^' ' 

where, for short, we have put 

a).(x) = Oi - X + fx - Ci - 2 j log ( 1 - ^. ) + ^ii^) + Hx - Oi). 

We next expand w,(x) in terms of powers of l/a,-. In the first place we 
have 



_a; + (x-a..-^)log(l-|) 



-X {^'-h') {^'^'-'-^^') 



+ 



Tlkix, Of) 



1 • 2oi 2 • 3a.2 kik + l)ai* ^ ai''+^ ' 

where r}{x, ai) is boimded. 

The function S,(x) is already a polynomial in l/o„ so that we need 
merely transform \(/{x — ai). 
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To find the desired expression for f{x — a.) let us formally expand the 
divergent infinite series 

Bi B, , Bz 



1 • 2(x - a.) 3 • 4(x - c.)' ^ 5 • 6(x - o.)^ ' 

as a series in IfOi. We obtain for the typical term the expression 

where the quantity in parentheses ends with (— lY'+^'^'^B^.+iyi, if s is odd, 
and in the term in x if s is even. We can simplify this result by making use 
of the Bemoulh polynomials. 
We know that 

<Pn{x)=x"-^x'^-' + (^l)B^x--'-(^l)B,x'^-^+ ..., 

the expression terminating with the term in a; or x^. The general term of 
the above expansion can thus be written 

1 f .^. s + 1 



sis + l)al 



I x'^^ - ^-^x- - ^.+i(x) + (- ly 5^1 



if we agree to put B(,+i)it equal to zero when s is even. 

Let us next choose the arbitrary integer t which occurs in \p{x — o<) so 
that 2t + 1 '^ k + 1. It follows then, by making use of the expression 
last obtained, that we can write 

t^(x - o.) = j-72^^ {x'' - X - <p2ix) -Bi}+ 273^2 1 x' - 2 a;" - <Pzix) J 

R(x,ai) 
■*■ a/+i ' 

where R{x, a.) is bounded. 

We now have all the quantities which occur in w,(x) written in temos of 
l/oi. Adding them we find the simple result 

where M{x, a.) = i?(x, a^) + R{x, at) is bovmded. 

The essential property of the expression which we have derived for Ui(x) 
is that the terms involving l/a.-, l/a.■^ • • •, l/cti*, do not contain the variable x, 
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and are, moreover, the first terms in the series which occurs in the asymptotic 
form of T{— Oi). This fact enables us to determine in a simple manner suit- 
able choices for pi, pi, • • • . 

Let us put 

Pi = r(- a.) = -\/2^(- ai)-'"~h(^''*\ 
where 

7? R (- 1) "* Sk+i , . 

^(,ai) - a< 1 . 2a, "^ 3 • Aa,^ " ' ^ k(k + l)a.* "*" a.*« ' 

e(a,) being bounded. 

With this selection the typical term of our infinite product is, for large 
enough i. 



e 



(- a,YT{- Oi) 

where M'{x, at) is bounded. 

What we have proved is this: As soon as the region in which x can vary 
is known we can find a positive integer I, and a positive constant K, such that, 
if Pi = r(— a.), the general term of our product is given by the last formula, 
in which \M'(x, a,)l < K, provided only that i > I. 

But we know that the series 

1 



Ci*+1 



is convergent, and we thus see at once that the selection we have made for 
Pi makes our product converge uniformly in any region where its factors 
are all analjrtic. 

We have then as a solution of the given equation the function 



*+> »n(^) 



This solution is seen to have no zeros in the finite plane, while the nature 
and situation of its singularities follow at once from the theory of the gamma 
fxmction. 

Let us take next the case where the zeros of E(x) are such that they Ue 
to the right of a line parallel to the axis of imaginaries. We obtain a 
solution by merely everjmrhere replacing the gamma function in the above 
result by the second gamma function; the second gamma function, T{u), 

being the solution of 

fix + 1) = xfix), 

that is asjmaptotic to Stirling's series in thejeft half plane; it is connected 
with r(w) by the relation T{u) = (1 - e"'"^-')r(M). 
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Consider now the case where E{x) is any entire function of finite genus. 
To apply the results already obtained we divide the zeros into two groups; 
those on or to the left of the imaginary axis forming one group, and those 
on the right the other. Let us call these two groups Ci, (h, • • •, and 61, bj, 
• • • , respectively. The solution is then found in a perfectly obvious way as 
the product of the solutions of two equations of the tjrpes we have considered, 

and is 

*+' *.w _ *+* *.(') 

Since the second gamma function, r(u), is zero for w = 1, 2, • • •, our 
solution is a meromorphic function, with its zeros at the points 

X = bi + n (i, n = 1, 2, ■ • •)» 

and its poles at the points 

X = tti — n (i = 1, 2, • • •, n = 0, 1, 2, • • •)• 

So far we have assumed explicitly that ^(x) is of finite order. Let us 
now suppose that this is no longer the case, but that the function g{x) is 
still a poljmomial. By making a modification of our reasoning exactly 
simil ar to that for finding the Weierstrass primary factors for the corre- 
sponding case we derive a solution, which is obtained from the one above, 
for the case of finite order, by making the summation which appears as 
the exponent of e in the ith factor include terms from n = 2 to n = i. 

We shall finally consider briefly the case where the coefficient of f{x) 
in (A) is a meromorphic function. The results which we have already 
obtained will enable us to show the existence of a solution provided a 
certain restriction is fulfilled. 

Suppose first that the coefficient is the reciprocal of an entire function. 
The substitution /(x) = l//(x) reduces the equation at once to the equation 
we have considered. This enables us to treat the general case; for we know 
that any meromorphic function is the quotient of two entire functions. 
Let us denote the extraneous exponential factor in these two functions 
by e^''-"'' and e"'^^', respectively. Then if gi{x) — g^ix) reduces to a poly- 
nomial we see that the solution can be written as the quotient of the solu- 
tions of two equations, each of which we know can be solved. The solution 
will be a meromorphic function, with its zeros and poles at known points. 

The class of entire and meromorphic functions for which we have obtained 
a simple and completely determined solution of (A) is seen to be a very wide 
one. 

Bloomington, Indiana, 
June, 1913. 



